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We propose a new means for quick and efficient acceleration of protons and/or electrons in rela- 
tivistic bulk flows. The maximum attainable particle energies are limited either by radiative losses or 
by the condition of confinement in the magnetic field. The new mechanism takes advantage of con- 
version of particles from the charged state (protons or electrons/positrons) into neutral state (neu- 
trons or photons) and back. In most cases, the conversion is photon-induced and requires presence 
of intense radiation fields, but under special circumstances the converter acceleration mechanism 
may operate via other charge-changing reactions, for example, inelastic nucleon-nucleon collisions. 
Like in the traditional, "stochastic" (or diffusive) acceleration models, the acceleration cycle in the 
proposed scenario consists of escape of particles from the relativistic flow followed by their return 
back after deflection from the ambient magnetic field. The difference is that the charge- changing 
reactions, which occur during the cycle, allow accelerated particles to increase their energies in 
each cycle by a factor much larger than 2 and usually roughly equal to the bulk Lorentz factor 
squared. The emerging spectra of accelerated particles can be very hard and their maximum energy 
in some cases is larger than in the standard mechanism. This significantly reduces the required 
energy budget of the sources of the highest-energy particles observed in cosmic rays, acceleration 
mechanism has a distinctive feature - it unavoidably creates neutral beams, consisting of photons, 
neutrinos or neutrons, whose beam pattern may be much broader than the inverse Lorentz factor 
of the relativistic flow. Also, the new mechanism may serve as an efficient means of transferring the 
energy of bulk motion to gamma-radiation and, if the accelerated particles are nucleons, inevitably 
produces high-energy neutrinos at relative efficiency approaching > 50%. 

PACS numbers: 98.70.Sa, 98.62.Js, 98.70.Rz 



I. INTRODUCTION 

The cosmic rays (CRs) have been studied for a long 
time (see, e.g., [l| for a review). There are many mod- 
els of their (likely non-uniform) sources, but explaining 
the origin of the highest-energy particles (> 10^° eV) 
observed in CRs is still a challenging problem 0. In 
general, two scenarios have been proposed for the par- 
ticle acceleration in astrophysical environments. One is 
the acceleration by electric field in geometries, where this 
field is not perpendicular to magnetic field-lines, for ex- 
ample, in the vicinity of magnetized rotating neutron 
stars 0. Very hard particle spectra may emerge in this 
way, but, because of the curvature losses inherent to such 
geometries, the upper limit to proton energy appears to 
be below 10^° eV. 

Another class of scenarios assumes gradual, "stochas- 
tic" (or diffusive) acceleration of charged particles 
through multiple reflections from inhomogeneities of 
magnetic field in the environments where large velocity 
gradients are present (see for a review). According to 
the generally accepted view (see, e.g., 0,0,0) this mech- 
anism works equally well for shocks and shear flows, and 
the emerging spectra of accelerated particles are such, 
that only a small fraction of the total energy budget is 
contained in the most energetic particles. 

Furthermore, there is a serious obstacle in achieving 
high energies of the accelerated particles - their diffusive 



escape. In the environments, where bulk velocities are 
sub-relativistic {v <^ c), it takes many reflections to in- 
crease the energy of accelerated particle twofold. In this 
case, the mechanism can work only in the largest ob- 
jects in the Universe, like galaxy clusters or radio lobes 
and knots in active galactic nuclei (AGNs), and only un- 
der optimistic assumption about the diffusion coefficient 
for particles in magnetic field. For example, Bohm dif- 
fusion leads to particle escape on the timescale R^/rgC, 
where i?o is the accelerator's size and rg the particle's gy- 
roradius, while the acceleration occurs on the timescale 
^ {c/v)'^rg/c. The maximum attainable energy is defined 
by equating these timescales: rg = {v/c)Ro <C i?o- It is 
of the order of 10^^ eV in the most favorable cases. 

Ultrarelativistic shocks and shear fiows are more 
promising, since the parameter u/c — > 1. Moreover, the 
energy gain per cycle, i.e., at each reflection from the 
shock or from the shear flow boundary, can approach 
the factor of ~ F^, where F is the bulk Lorentz fac- 
tor (measured in the upstream fluid frame in the case 
of shock). Such a scenario was suggested, e.g., as the 
dominant acceleration mechanism in gamma-ray bursts 
(GRBs) Illlli. However, in the standard acceleration 
theory, the energy gain of F^ occurs only at the flrst 
cycle, while the subsequent ones result in energy gain of 
~ 2 each US [ill. 

The reason is that the shock catches up to the reflected 
particle (or the particle crosses the boundary of the shear 
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FIG. 1: The acceleration cycle in the standard mechanism 
for a shock (a) and for a shear flow (b). The thick solid line 
shows the particle's trajectory. The magnetic field is per- 
pendicular to the picture plane. The locations of the shock 
at the moments of particle escape from the shock and sub- 
sequent catch-up are shown as double lines. The shear flow 
boundary is shown by thin dotted line. 



flow) when its trajectory makes an angle of ~ 2/r to the 
shock normal (or bulk velocity vector) , while keeping the 
energy gain would require isotropisation. 

The smallness of this angle is apparent for the case of 
shear flow (see Fig. lb). For shocks it is a consequence of 
relativistic motion (Fig. la). Indeed, the component of 
particle velocity along the shock normal becomes smaller 
than the velocity of the shock itself as soon as the an- 
gle between the particle's momentum and the shock nor- 
mal grows larger than l/F. There is no apparent way to 
isotropize reflected particles unless a special structure in 
the magnetic field ahead of the shock is introduced, e.g., 
a counter-propagating shock. 

However, a means to circumvent these limitations ex- 
ists ~ it can be done by switching particle's charge 
on and off at right times. Paradoxically, interactions 
with photons, which have been always treated as dissipa- 
tive processs leading to degradation of particle energy, in 
fact play positive role: they allow (through the charge- 
changing particle conversion) to keep the energy gain 
up to the largest particle energies. There are also other 
types of conversion reactions, which we briefly discuss 
in the following section. One of the suggested below 
charge-changing schemes (electron-photon-electron reac- 
tion chain in GRB shocks) is independently considered 
in [3. 

In order to outline the general picture, we intentionally 
skip some details, which are not essential for the proposed 
acceleration mechanism, but may change its quantitative 
characteristics. In particular, we assume that the mag- 
netic fleld is either chaotic (turbulent) or uniform with 



field lines perpendicular to both the momentum of accel- 
erated particle and the velocity of the flow, and treat 
shocks and shear flows as one-dimensional discontinu- 
ities. 

We use the convention F^, oc i''? to define the spectral 
index q, where F^, is the energy flux per unit frequency 
interval. 



II. THE CONVERTER ACCELERATION 
MECHANISM 

Two basic types of photon-induced conversion involve 
nucleons or electrons/positrons. Both cycles consist of 
two reactions: 



p + J ^ n + TT 
for proton cycle, 

-I- 7 — > + 7' 



and 



and 



n + ■y ^ p + TT 



7' + 7 ^ 6+ + e 



(1) 



(2) 



for electron cycle. Here p, n, n^, and denote proton, 
neutron, charged pions (positive and negative), positron 
and electron, respectively; 7 is a low-energy background 
photon and 7' the high-energy comptonizcd photon. 

The second of reactions from the electron cycle © 
has a kinematic threshold — 2 rrieC^ in the center-of- 
momentum frame, where nie is the electron mass. Ef- 
fectively, the first reaction also has the threshold ~ Ag, 
since at lower energies of incident photons the efficiency 
of energy transfer to the comptonizcd photon becomes 
much less than unity. The reactions proceed differently 
depending on the background photon spectrum. Soft 
spectrum blocks the electron cycle as the fraction of en- 
ergy transferred to a comptonizcd photon is too small 
and there are few target photons sufficiently energetic 
for the second reaction of the electron cycle. An example 
of soft spectrum could be a power-law with spectral index 
g < — 1 or a narrow-band spectrum, like black-body or 
line emission, with typical photon energy e <^ m1c^/ee, 
where £e is the electron energy. For intermediate spectra 
(e.g., power- laws with indices — 1 < g < 1 or narrow- 
band spectrum with e ^ m'lc'^/ee), the comptonizcd pho- 
ton takes about 1/2 of the electron (positron) energy, 
and in the consequent pair production event this energy 
is divided nearly in equal parts between the daughter 
electron and positron. The cross-section in both pro- 
cesses is (Te,-Y ~ 10^^^ cm^. For hard spectra (g > 1 
or e ^ TOgC^/£e), these reactions proceed in the deep 
Klein-Nishina regime, i.e., the comptonization and pair- 
production cross-sections decrease inversely proportional 
to the square of the center-of-momentum energy (their 
ratio is 1:2) and almost all the energy of interacting par- 
ticles is transferred to one of the daughter particles. In 
effect, the energy losses for the combined electron/photon 
particle become very gradual. This case is the closest 
to the pure conversion (the probability of charge change 
Pc = I) without accompanying energy losses provided the 
synchrotron emission is negligible. 
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The reactions from the proton cycle 1^ (see, e.g., 
have kinematic threshold of (m^c^ + mJc^/2 mTv) — 
150 MeV in the nucleon rest frame (m7r ~ 140 MeV/c^ 
is the charged pion mass and rriN ~ 940 MeV/c^ the 
nucleon mass). Side by side with the reactions 1^ pro- 
ceed other photopionic reactions with formation of neu- 
tral pion, which preserve nucleon's charge. They have 
roughly the same cross-section and should be considered 
as a background energy losses. The total photopionic 
cross-section rapidly increases with energy of incident 
photon and reaches maximum of ct^ ~ 6 x 10^^* cm^ 
at Ap ~ 340 MeV, which corresponds to formation of A- 
resonance and should be considered as effective thresh- 
old. Well above the resonance energy, the cross-section 
decreases and levels off at ~ 10~^^ cm^. The probability 
of charge change in a photopionic reaction is Pc — 1/3 at 
the resonance and Pc — 1/2 at the plateau. The inelas- 
ticity is ~ 0.2 and about 0.5, respectively. 

A competing photon-induced reaction is the process of 
creation of an electron-positron pair by a photon interact- 
ing with the electric field of a proton 7 ^ p + e~ +e~^, 
which has the cross-section ~ 5 x 10~^^ cm^ and inelas- 
ticity ~ 10"'^. With the decrease of spectral index q of 
target photon field, this process becomes an increasingly 
important energy loss channel, and at g ~ —1.55 the dif- 
ference with the photopionic processes in the inelasticity 
and cross-section is exactly balanced by larger number of 
target photons (thanks to lower threshold) . Anyway, the 
p-l-7— >p + e^+e+ process can be neglected for spectra 
with q > —1.5. 

In dense environments and at relatively low nucleon 
energies, i.e., in the case where there are few target pho- 
tons, the proton cycle proceeds through inelastic nucleon- 
nucleon collisions, for example, 



p + p n + p- 



and 



n+p p + p + n 



(3) 



The kinematic threshold for these reactions is rriT^c^ in 
the center-of-momentum frame, and the cross-section at 
energies ^ m^rC^ is ~ 3x 10^^® cm^. The acceleration via 
inelastic nucleon-nucleon collisions could be important in 
GRB internal shocks, where the required column density 
of ~ 10 g/cm^ is achieved [l5l |. 

In both proton and electron cycles, one can consider 
an accelerated nucleon or electron/positron as a particle, 
which has both charged and neutral states. The accel- 
eration cycle consists of three steps (see Fig. 2). First, 
a charged particle in relativistic flow is converted into 
neutral state (point 1). Then, experiencing no influence 
from the magnetic field, it freely leaves the flow and prop- 
agates into ambient medium much further than if it were 
charged. Second, a transition from neutral to charged 
state occurs (point 2), which may be the spontaneous 
neutron decay. 

At this moment, the particles in the laboratory frame 
preserve beaming with the opening angle of 1 /F, which 
they had in the neutral state. The initial handicap al- 
lows particles to be deflected by an angle ^ 2/F be- 
fore the encounter with the relativistic flow. The angular 
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FIG. 2: The acceleration cycle in the converter mechanism for 
a shock (a) and for a shear flow (b). The particle's trajectory 
is shown by thick dotted line (neutral state) and thick solid 
line (charged state). The magnetic field is perpendicular to 
the picture plane. Numbered are the moments of particle con- 
version into neutral state, transition from neutral to charged 
state, and subsequent return to the flow. The locations of 
the shock at the corresponding moments are shown by double 
lines. The shear flow boundary is shown by thin dotted line. 



spreading of the particle beam in the laboratory frame 
means energy gain in the flow comoving frame, which is 
much larger than 2 and amounts to F^ in the case of full 
isotropisation. At the third step, the particles return to 
the flow (point 3) and their isotropisation in the comov- 
ing frame translates into the resulting energy gain in the 
laboratory frame. 

The main parameter, characterizing the efficiency of 
the converter mechanism, is the optical depth for inter- 
actions of accelerated particles (protons/neutrons or elec- 
trons/photons) r — anD, where n is the number density 
of target particles (photons or nucleons) and D the ac- 
celerator's size, both measured in the comoving frame, a 
the relevant cross-section. 

The optical depth is geometry-dependent. In the case 
of a continuous outflow or a shock, produced by a cen- 
tral engine and subtending an angle > 1/F, one has for 
photon-induced reactions 



4 nRce^ 



«2 



(4) 



Here L is the apparent luminosity per logarithmic fre- 
quency interval at photon energy = 2mc^A/(eO^), 
where the interaction with target photons is the most 
efficient, £ the energy of accelerated particle, R the dis- 
tance from the central engine, A and m are the threshold 
and the mass of the particle for one of the possible cycles. 

The beaming angle of target radiation field is O ~ 1/F 
in the case where it is produced within the jet or by the 
shocked gas, and 0^1 for the emission from broad-line 
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regions in AGNs and the radiation scattered in the inter- 
stehar medium around GRBs. Intermediate cases, where 
1/r < < 1, are also possible. They include, for ex- 
ample, the radiation from inner parts of accretion discs 
in AGNs or the emission which accompanies the acceler- 
ation process (as explained in Sect. V). If the source of 
target radiation is transient, like a GRB, with duration 
less than RQ^/c, then its luminosity should be replaced 
by Ec/{R@'^), where jLdt. 

The most favorable for the converter mechanism con- 
ditions exist in AGNs and GRBs, where both ultrarela- 
tivistic flows and intense radiation fields are present (see 
[TgI and for reviews). Other objects with relativis- 
tic outflows, e.g., stellar-mass microquasars, may also 
be able to accelerate particles via the converter mech- 
anism. The only two prerequisites are sufficiently high 
conversion probability (which we specify below) and bulk 
Lorentz factor F 3> 2 in order to compensate the energy 
losses caused by conversion. 

Because of relatively large cross-section, the optical 
depth for photon-induced reactions is not a bottleneck 
for the electron cycle, but it could be a limiting factor 
for the proton cycle. Let us estimate the optical depth 
for photopionic reactions in three cases. For AGN inner 
jets we obtain (taking into account only comoving photon 
fields with 6 = 1/F) 



VlO^^erg/sy VlQiSeV/ 



10 V /lO^^cm 



(5) 

where the apparent luminosity per logarithmic frequency 
interval L(£*) depends on e. The radiation from AGN 
broad-line regions creates the optical depth 



T - 5 X 10-2 



L 



10^4 erg/s 



10eV\ /lO^^cm 



(6) 



which does not depend on the particle's energy: for all 
particles with energy e > 2mc^A/e ~ 5 x 10"'^^ eV the 
number of target photons is essentially constant because 
of their relatively narrow spectral distribution. In GRBs, 
the optical depth due to comoving photons is 



T ~ 3x10" 



IQ52 gj-g 



100 \ / lO^^cm 



R 



(7) 

So, the conversion probability Pcn = [1 ~ exp(— pc^)] in 
the proton cycle is usually, though not always, much 
smaller than unity. 

It should be noted, that effectively the probability of 
conversion of a neutron into a proton always exceeds 

^(..) ^ ^ ^ 3 X 10- C-^) (^) (8) 

because of the neutron decay. Here i„ ~ 900 s is the 
lifetime of free neutron. The spontaneous decay of free 
neutrons is important at small energies, especially during 



the first acceleration cycle, while at large energies the 
photon-induced conversion is more efficient. 

The expressions (0 and Q cover all physically 
different situations. For example, one may use Eq. © 
to estimate the optical depth in microquasars, where the 
target photons are produced by a hot corona having the 
size R ^ 10* cm and the luminosity L ^ 10"^^ erg/s. The 
result is T ~ 0.1 at target-photon energy of ~ 1 keV, 
so that the proton cycle may operate in microquasars at 
energies > 3 x lO^'* eV. It is possible, therefore, that some 
contribution to the galactic CRs around the knee comes 
from microquasars. In any case, we see no problem in 
realization of the electron cycle in such objects. 



III. ENERGY GAIN 

Assuming that momenta of particles are isotropized 
in the comoving frame upon their encounter with the 
relativistic flow, one gets the average energy gain per 
cycle 



9 ■ 



(9) 



where 9 is the angle between the particle's momentum 
and the flow's velocity (deflection angle) at the moment 
of encounter. If the deflection angle is small, 9 <^ 1, then 
it grows linearly with distance i (travelled by the particle 
in the charged state) 



(10) 



in uniform magnetic field, while in chaotic field it behaves 
as 



(0) 



(11) 



on average. Here 9q is the initial deflection angle, mea- 
sured at the time of conversion from neutral to charged 
state, £c < i the turbulence scale of the magnetic field, 
and Vg the gyroradius of the particle, calculated as if the 
field were uniform. Note, that the turbulence scale is im- 
plicitly defined by Eq. Hll|l and, therefore, may depend 
on the particle's energy. 

In the case of acceleration at the shock front, the 
shock catches up to the particle when the displacement 
of the particle along the shock normal, after it crossed 
the shock, becomes equal to the distance traveled by the 
shock front, i.e., 



(4 +^) = 4cos6io + / cos9d£', (12) 



where £o is the distance travelled by the particle in the 
neutral state after it left the relativistic flow, and we 
neglected the difference between the particle's velocity 
and the velocity of light. 
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In the uniform magnetic field, the deflection angle is 




provided rg/T <C ^ ^g^^- For smaller initial dis- 
placement, one arrives to the result of standard theory 
~ 1/r, whereas for larger displacement the deflection 
angle is 6* ^ 1. In the rightmost part of Eq. ^1?>\ we as- 
sume that [l-T'^el) - 1 to simphfy the algebra. This 
corresponds to ignoring particles, which escape propagat- 
ing (in the comoving frame) nearly parallel to the shock 
plane. They arc neither numerous nor energetically im- 
portant. 

If one substitutes by R (or by YD in the case of a 
small blob with the size D < R/T, ejected by a central 
engine), then the applicability limits give two critical en- 
ergies (both measured at the end of acceleration cycle - 
point 3 in Fig. 2): 

£i = TeBR and £2 = eBR, (14) 

where e is the charge of the accelerated particle and B 
the magnetic field strength. The acceleration proceeds 
with the maximum energy gain of ~ up to the energy 
62, whereas above the energy ei there is no advantage 
over the standard mechanism in the energy gain. At the 
same time, Si is equal to the maximum energy, achievable 
in the standard mechanism. 
In the chaotic magnetic field 




for Vg/T^ <C tc^o F^r-g, and the critical energies are 
equal, 

£1 = £2 = FeSy!R4- (16) 
In the case of shear flow, the catch-up condition reads 

4 sin 6*0 = / sin 61 d^'. (17) 
Jo 



The deflection angle is 




for chaotic field (r^/F^ < < r^T). 

There is a subtlety in realization of the converter 
mechanism if acceleration takes place at a shear flow 
boundary. If the charged particle reappears at a dis- 
tance greater than Vg (uniform magnetic fleld) or r'^/ic 



(chaotic field) from the fiow boundary, then it should 
drift or diffuse back. It takes a time of the order of 
td — R"^ / (Trgc), in the least favourable case of quasi- 
uniform magnetic field with characteristic spatial scale 
of ^ R, or td — R'^ic/ {T'^r'^c), in chaotic magnetic field. 
In both cases we assumed ^gsin^o = R/^- The solution 
of diffusion problem with a sink (the shear flow) in the 
half-space shows, that all particles eventually return to 
the flow, provided it persists for sufficiently long time. 
In reality, some of them are lost because the shear fiow 
does not occupy the entire half-space, but the losses are 
negligible if the spatial extent of the shear flow boundary 
is much larger than ^osin^o- This condition, for exam- 
ple, is satisfied for a conical jet with the opening angle 
> 1/F. 

The synchrotron losses during the time td might be a 
more serious problem. One can neglect them if the energy 
loss rate £ = 4/9 {e^ / mc^Y B'^ {e / mc^Y c multiplied by td 
is less than the energy of accelerated particle. In the case 
of quasi-uniform magnetic field, this condition turns out 
to be energy-independent and, after simple algebra, one 
gets 

R»-a^.—2{^)'^ i-^-' (20) 

Here i?opt is the optimal size of an electromagnetic ac- 
celerator 0], which is defined to minimize the amount 
of energy contained in the electric and magnetic fields. 
Let us recall, however, that the optimal size is originally 
defined for the accelerator itself (in this case ~ for the 
interior of the shear flow), while Eq. (|20|l has to do with 
the external magnetic fleld. 

In a turbulent magnetic field, the synchrotron losses 
arc the most significant for the least energetic particles, 
i.e., those having rg ^ £c- The losses can be neglected if 



where the optimal size is the function of critical energy 
£2 , which is defined as if the magnetic field were uniform. 

In order to be in agreement with the observed diffuse 
gamma-ray background, a typical source of the highest- 
energy CRs must satisfy the condition R > 0.2i?opt Hi; 
so that the synchrotron losses during diffusion do not lead 
to any further resrictions. The aforementioned statisti- 
cal arguments are not applicable to individual abberrant 
accelerators, for which Eq. H20|l or Eq. (|21|l should be 
considered as an additional limit. 

The critical energies in the case of acceleration by a 
shear flow are 

£1 = £2 = TeBR (22) 

for uniform magnetic field and 

Si^TeB^/Rfc and £2 = F^e^V^ (23) 

for chaotic field. It is interesting to note, that in the lat- 
ter case £2 > £1, i.e., the maximum energy achievable in 
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the converter mechanism appears to be T'^/'^ times larger 
than that in the standard mechanism. However, a par- 
ticle can only attain this energy if it enters acceleration 
cycle with the energy e ~ £2/r^- A particle, entering the 
cycle with larger energy, will end up with smaller energy, 
in contrast to the other three cases. 

Of course, every acceleration mechanism must obey 
the fundamental restrictions and constraints of classical 
electrodynamics 0. Among them is the Hillas criterion 
[TsL [l9l | generalized for relativistic bulk flows, which is 
simply equivalent to the condition £ < ei, where ei is 
the first critical energy from Eqs. (|14|l and (|22|1 . This 
is not a chance coincidence, but rather the consequence 
of the implied - uniform - magnetic field configuration. 
In a turbulent magnetic field, the maximum attainable 
energy is always lower than the fundamental limits, but 
the converter mechanism is less affected by the turbu- 
lence. Indeed, the magnetic field should be considered 
chaotic if its turbulence scale is less than the gyroradius 
of a particle at the beginning of the acceleration cycle. It 
means £c < ^g(£i) for the standard mechanism, but only 
< ?'g(£2/r^) for the converter mechanism. Thus, a 
particle accelerated in the standard way feels turbulence 
starting from larger scale, and a situation is possible, 
where the magnetic field should be treated as chaotic in 
the standard mechanism while being essentially uniform 
for the converter mechanism. The inequality £2 > £1 
from Eq. I|23|l is, in fact, the consequence of slower de- 
crease of critical energy with decreasing £c for the con- 
verter mechanism. 

So far, the analysis was limited to the case of negli- 
gible synchrotron losses and r <^ 1. The larger is the 
probability of conversion, the easier is acceleration in the 
converter mechanism. But if t approaches unity, then 
the critical energies decrease: they are still given by Eqs. 
((Tljl . (|22J), and (ESJl, but the value of R should be 

replaced by the (energy-dependent) mean free path of the 
accelerated particle. Analogously, the synchrotron losses 
also decrease the critical energies; in this case R should 
be replaced by the radiation length. 

In both cases the energy share of the accompanying 
radiation is non- negligible or even dominant, although 
the status of conversion losses is qualitatively different 
for the proton and electron cycles (see Sect. V). 
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FIG. 3: The particle distribution (solid line), resulting from a 
monoenergetic injection. The dips in the distribution preserve 
themselves until the width of the injection spectrum is larger 
than logF^ in logarithmic units. 



where the particle's energy is simply proportional to the 
energy gain, e = gso. It is a power-law distribution, but 
all the particles with displacement larger than the ap- 
plicability hmits of Eqs. and have 
the constant energy gain of ~ F^, so that in general a 
delta-function is added to the resulting particle distri- 
bution at its high-energy end. In their power-law parts, 
the distributions emerging in uniform and chaotic field, 
respectively, are 



dN 



OC £ 



1/2 



and 



dN 



de de 
for acceleration at the shock front, and 



OC £ 



dN 
de 



const 



and 



dN 
de 



OC e 



1/2 



(25) 



(26) 



for acceleration at the boundary of shear fiow. These 
spectra are extremely hard. In practice, they can be con- 
sidered delta-functions, so that the resulting distribution 
is defined mainly by the spectrum of injected particles. 

After a few cycles, a seesaw-shaped spectrum (Fig. 3) is 
formed from the initial monoenergetic distribution. The 
number of particles emerging from a cycle, N' , is related 
to the number of particles entering it as N' = kN, where 



IV. RESULTING PARTICLE DISTRIBUTION 

Let us consider a monoenergetic {e = £0) beam of par- 
ticles in the neutral state, escaping from the relativistic 
fiow. If the probability of conversion per unit length. A, 
is constant (i.e., we neglect the exponential decrease in 
the number of particles, assuming r <C 1), then the dis- 
tribution of charged particles over initial displacement 
dN/d£o = A translates into 



dN dN dio , f„ dO 
— — OC — — — OC A o— — 

de d£o dg \ dto 



(24) 



■ _l J_ Pen 



(27) 



is the overall probability of conversion, which normally 
occurs twice per cycle (points 1 and 2 in Fig. 2), mul- 
tiplied by the probability of escape from the relativis- 
tic fiow after the first conversion Pcsc, which is the rel- 
ative number of particles moving upstream with respect 
to the shock front or the shear fiow boundary. Assum- 
ing that the particle distribution is isotropic just before 
the conversion into neutral state, one obtains Posc = 1/3 
for a strong relativistic shock, where the shocked fiuid 
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moves at the speed c/3 away from the shock plane, and 
Pesc = 1/2 for a shear flow, where the fluid velocity is par- 
allel to the boundary. Within one acceleration cycle, the 
probabilities p"cn and p"cn are of the same order, except 
for the case where pin is limited by Eq. 

There is a simple case, where k remains constant in 
the course of acceleration. This happens in both electron 
and proton cycles if the spectrum of target photons is a 
power-law with spectral index g = or (for the proton cy- 
cle only) in the case of a narrow-band target-photon field. 
Then the envelope of the particle distribution shown in 
Fig. 3 is a power-law dN/de oc £~", and the index a can 
be obtained in the following way. 

Since the particle energies are related as e' = ge, we 
find that 

where primed is the distribution emerging from the cycle. 
Taking logarithm of both parts, we obtain 

, ^ k\ Ink , ^ 

-aln5 = ln - ^ a = l-- — . (29) 

If the probability of passing through the acceleration cy- 
cle k is larger than the inverse energy gain 1/g, then the 
spectral index is a < 2, i.e., most of the energy content is 
at the high-energy end of the distribution. Both AGNs 
and GRBs can satisfy this requirement (see Eqs. [SJ 
and [7|l and, therefore, can be efficient producers of the 
highest-energy cosmic rays. 

In the energy range, where the efficiency condition for 
the converter mechanism is satisfied, this mechanism is 
the dominant source of accelerated particles - just be- 
cause it provides a spectrum harder than the one result- 
ing from the standard mechanism. At the same time, 
there is no actual threshold in the conversion probability: 
the converter mechanism can function even at Pcn 0, 
but the number of the accelerated particles in this case is 
extremely depleted at high energies. Quantitative studies 
of the emerging CR spectrum must take into account the 
dependence of the conversion probability on the energy of 
accelerated particles. Hence, the precise solution can be 
obtained only in the self-consistent approach, which in- 
cludes the effect of the accompanying emission, discussed 
in the following section. 



V. THE ACCOMPANYING EMISSION 

Acceleration of particles via the converter mechanism 
is inevitably accompanied by gamma-ray and - in the 
case of nucleon conversion - by neutrino emission. The 
conversion losses in the electron cycle, which are the re- 
sult of non-optimal inverse Compton scatterings, do not 
lead to irreversible energy drain. Instead, all the energy 
remains in photons or electrons (positrons) and, as they 



both can in principle participate in the acceleration cy- 
cle, may be used to inject new particles for acceleration. 
Thus, even in the case r ^ 1 the influence of the conver- 
sion losses is reduced to the decrease of maximum attain- 
able energy and the increase in the level of synchrotron 
emission by rising the number of particles involved in the 
acceleration. 

On the contrary, the conversion losses in the proton 
cycle act as a true energy sink: the energy spent for the 
pion production never comes back to the accelerated nu- 
cleons. Neutrinos, e~e"'"-pairs and gamma-rays, which 
are the decay products of pions, are copiously produced 
as byproducts of the proton acceleration cycle. Neutri- 
nos carry away about one-half of the energy of acceler- 
ated nucleons and, because they freely escape from the 
acceleration site, their spectrum copies that of the nu- 
cleons, but scaled down in energies by about an order 
of magnitude. The gamma-rays and pairs, which carry 
approximately the same energy as neutrinos, are repro- 
cessed through the electromagnetic cascade: the photons 
are absorbed in two-photon pair production process, elec- 
trons and positrons cool in the magnetic field, producing 
synchrotron radiation, i.e., another generation of pho- 
tons, etc. 

The hard spectrum of primary electrons and positrons 
means that they form a standard cooling distribution, 
dNe/dje oc 7~^, and their synchrotron spectrum is a 
power-law with (photon) index —3/2. The synchrotron 
photons spawn another generation of pairs, which cool to 
form a distribution with index —5/2. The cascade com- 
prises several steps like this, so that the photon spectrum 
eventually converges to the power-law with index —2. 
However, at low energies the acceleration site is transpar- 
ent for the photons, and therefore pairs are not injected 
below some energy, which leads to a break in the spec- 
trum; below the break the photon spectrum preserves its 
original spectral index —3/2. The location of the break 
can be found in a self-consistent way: 

...-0.5-?^(-ii^)\ (30) 

1 TTleC \ 2 mgC^ / 

where h is the Planck constant. The threshold photon 
energy £tr is defined so that the optical depth for two- 
photon pair production, 

"-^^*') = ^"-2Tw(,^j ' 

is equal to unity, where a-y-y is the cross-section of two- 
photon pair production and the optical depth for pho- 
topionic reactions at photon energy e*. Straightforward 
calculations yield 

e,.^0.5^J^^^f^. (32) 

Substituting B = 10^ G, En = 10^^ eV (the energy 
of accelerated nucleons), F = 1000, and r^r = 10~^, - 
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parameters, reasonable for GRBs, - we obtain the break 
energy etr ^ 100 keV, which is similar to the really ob- 
served in GRBs spectral features. Within this picture, 
however, there is no simple way to explain the observed 
relative stability of £br whereas the parameters entering 
Eq. (|32|l . especially Ttt, may vary by orders of magnitude. 

Unlike the conversion losses, the synchrotron emission 
is not tightly related to the acceleration process as such. 
It may have a negligible effect, especially for protons, but 
also may be the main energy loss channel. The detailed 
analysis of the properties of accompanying synchrotron 
radiation is beyond the scope of this paper, but we point 
out two distinctive features. 

First, the maximum energy of synchrotron photons for 
the converter mechanism is times larger than for the 
standard one. The existence of such an energy limit is 
easy to see for the standard mechanism (following the 
arguments of '53]). The acceleration cycle in this case 
lasts ~ ''ff/c and the energy increment is ^ £, which 
gives the acceleration rate i ~ £c/rg. The maximum 
admissible rate of synchrotron losses is just the same, 
so that the particle's energy is limited by the following 
inequality: 




Then, simple calculation yields the maximum energy of 
synchrotron photons (ignoring relativistic dipole radia- 
tion caused by small-scale inliomogeneities of magnetic 
field): 

e('^--)^0.5—(^Y ^ ^mc' ~ 137 mc\ (34) 

Doppler boosting gives additional factor F, and in a 
turbulent magnetic field with the spatial scale of the 
turbulence less than — m^c^/ve^B^ another factor 
{£c/ecT^^ apphes. 

The same reasoning is valid for the converter mech- 
anism, in which the cycle duration is > rg/c and the 
energy increment is < F^e. Consequently, the analog of 
Eq. H33|l gives F^ times larger limit for the energy of ac- 
celerated particle, which translates into factor F* in the 
expression for the energy of synchrotron photons. When 
the accelerated particle enters the relativistic flow be- 
ing close to the limiting energy, the synchrotron emission 
is so efficient, that the particle loses almost all its en- 
ergy before it is deflected by an angle ^ 1/F. Thus, 
the resulting synchrotron emission is beamed backwards 
in the flow comoving frame. In the laboratory frame it 
appears redshifted by the factor F, in contrast to the 
standard mechanism, in which the synchrotron emission 
is blueshifted in the laboratory frame by the same fac- 
tor F. Thus, for an observer resting in the laboratory 
frame, the maximum energy of synchrotron photons ac- 
companying the converter-acceleration is by the factor of 
~ F^ larger (recall the difference of ~ F^ in the comoving 
frame), than in the standard mechanism. Moreover, this 



highest-energy synchrotron radiation is quasi-isotropic in 
the laboratory frame, which is another distinctive feature 
of the converter mechanism. The latter phenomenon has 
a nature similar to the effect of beam-pattern broadening 
for the inverse Compton radiation of electrons in front of 
relativistic shock [21|. Generally speaking, the converter 
mechanism makes neutral beams of all kinds (photon, 
neutrino and neutron beams) broader than 1/F, so that 
they can be seen even if the jet that produced them is 
not pointing towards the observer. 

The accompanying electromagnetic emission of any 
origin can itself provide photons for the conversion reac- 
tions, which may give rise to radiative instabilities analo- 
gous the instability driven byp+7 p+e^ process, 
discussed in [2^ . 



VI. DISCUSSION 

The key parameter of the converter mechanism - the 
probability of conversion - varies from one cycle to an- 
other or even within cycles. It increases in the course of 
acceleration (because a more energetic accelerated parti- 
cle interacts with less energetic and more abundant target 
particles) and with increase of the angle between the par- 
ticle's trajectory and the velocity of the flow. The flrst 
effect is negligible for nucleon-nucleon collisions, while 
the second is partially or even completely compensated 
by smaller distance travelled at large angles. There is, 
however, a general trend of increase of the probability of 
conversion towards higher energies of accelerated parti- 
cles; this may even block further acceleration if pcn ap- 
proaches unity. Nevertheless, with the expansion of the 
flow the density of target particles drops, leading to the 
decrease of the probability of conversion, and the accel- 
eration resumes. 

It means, that the converter mechanism is capable of 
self-tuning. The only thing which is required for this 
mechanism to be efficient is that the probability of pass- 
ing through the acceleration cycle k is larger than 1/F^ 
somewhere along the flow. It does not matter how large 
is the probability: in the case A: — > 1 the particles are 
preaccelerated in the region of high optical thickness and 
then are further accelerated when (or where) the flow be- 
comes optically thin for them. Of course, the question 
how much energy is wasted for conversion losses during 
this intermediate phase is open. 

On the other hand, the condition k > 1/F^ places 
some restrictions on the sources where the proton cy- 
cle is readily realizable starting from thermal protons. 
If a source is known to be bright at hard gamma-rays, 
one can conclude that the source's compactness, calcu- 
lated for the corresponding target photons, is low, so 
that the probability of proton-to-neutron conversion is 
k ^ cTTr/aj ~ 3 X 10~'^. Thus, the source should have 
F > 20 - a condition, satisfied by GRBs and many AGNs. 
Moreover, this condition is not obligatory for the sources, 
where the converter mechanism starts from preacceler- 
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ated high-energy protons, which require only low-energy 
target photons. 

We expect that the converter mechanism never oper- 
ates alone and the standard Fermi-type mechanism com- 
petes with it. Because the particles spend some time in 
the neutral state, the converter mechanism has smaller 
acceleration rate at low energies. But, close to the lim- 
iting energy (ei), the durations of acceleration cycle in 
both mechanisms approach the same value i?/c, so that 
the average acceleration rates are roughly equal. Except 
for the case of acceleration in a shear flow with chaotic 
ambient magnetic field, the mechanisms have just the 
same absolute, i.e., ignoring radiative losses, energy limit 
for the accelerated particles. However, in the converter 
mechanism, this energy is attained in much fewer steps, 
with potentially much more particles survived. More- 
over, the particles - when converted into neutral state 
- can escape from regions located deep inside the rela- 
tivistic flow, which further reduces irreversible particle 
losses in the downstream. So, at the highest energy part 
of the distribution almost all particles are produced by 
the converter mechanism, regardless of its performance 
at low energies. 

Now let us consider the phase, when produced CRs 
leave the accelerator. Within the framework of the con- 
verter mechanism, the particles can escape in the form 
of neutral beam - an easy way, which causes no problem 
and requires nothing but a sufficiently high conversion 
probability. In the standard scheme, the escaping par- 
ticles are charged and inevitably must form an expand- 
ing turbulent outflow. This causes considerable adiabatic 
losses, which can hardly be controlled. Another advan- 
tage of the converter mechanism is its greater tolerance 
for non-uniform magnetic fields, as discussed in Sect. III. 
Since the magnetic field turbulence has strictly negative 
effect on the maximum attainable energy, the converter 
acceleration mechanism may have a larger cut-off energy. 

Unlike the standard mechanism, the converter mech- 
anism in many cases does not need any special particle 
injection. By definition of acceleration cycle, even the 
particles resting either in the laboratory or in the flow 
frame serve as injection as soon as they are converted into 
neutral state. Sometimes however, e.g., in AGNs, the 
probability of conversion is too small at low energies to 
provide sufficient injection. Then, the converter mecha- 
nism comes to depend on the standard mechanism, which 
produces preaccelerated particles (at energies ^ 10^^ eV 
in the case of AGNs). 

The photopionic processes in AGNs and GRBs lead- 
ing to copious production of neutral particles have been 
extensively discussed in the literature (e.g., 0,11^), in 
particular in the context of predictions of fluxes of high 
energy neutrinos from these objects. The main reason- 
ing was that these objects are postulated to be efficient 
sources of high-energy CRs and they are surrounded by 
dense photon fields. Therefore, nucleon-photon interac- 
tions may lead to detectable neutrino fluxes. Our point 
of view is essentially different: if one wants to accelerate 



CRs in these objects to highest energies, then the con- 
verter mechanism is likely the most efficient way to do 
this, and the neutrino emission appears as a natural and 
unavoidable byproduct of this mechanism. 

Further questions, absent in the test particle problem, 
arise in the self-consistent approach. In particular, the 
converter mechanism strongly alters hydrodynamics of 
relativistic shocks or shear flows. We have shown that 
they can accelerate nucleons with the efficiency of energy 
transfer approaching 100%. At such a high efficiency, 
almost all the available energy is transferred to the end 
of the particle distribution, and hence the ultra-energetic 
particles contribute most of the inertia of the relativistic 
flow (shocked gas). Under such circumstances, the idea 
of a shock as a discontinuity becomes meaningless, since 
the gyroradius of dynamically most important particles 
is comparable to the size of the system. 

This limits applicability of the test particle approach, 
since it is only valid for the particles whose gyroradius 
is larger than the width of the shock or the shear flow 
boundary. The similar problem exists also in the stan- 
dard mechanism, which is, moreover, more sensitive to 
it. Indeed, in the converter mechanism an accelerated 
particle may cross the shock or the shear flow bound- 
ary while being in the neutral state, hence relaxing the 
requirements to their sharpness. 

The absence of a true shock is an obstacle for the elec- 
tron cycle, which would otherwise be a very efficient way 
of direct acceleration of electrons. When the mean free 
paths of photons and pairs are small compared to the 
spatial scale of velocity gradient, the direct acceleration 
of electrons probably cannot rival the efhciency of pro- 
duction of secondary pairs via nucleon acceleration. One 
of the implications is that at least in some, particularly 
bright, GRBs the main source of the observed gamma-ray 
radiation is the pion-induced cascade. 

Energetic neutrons (or photons) escape from the rela- 
tivistic flow and their decay or interaction products can 
disturb the ambient medium to the extent that it starts 
to move with an ultrarelativistic speed. An apparent con- 
sequence of this is the distortion of beam pattern of the 
source, while the influence on the limiting energy of the 
accelerated particles needs further investigation. On one 
hand, the decrease in the Lorentz factor contrast lowers 
the limiting energy. On the other hand, enhancement of 
the ambient magnetic field, either via compression in the 
shock or by means of various instabilities triggered by 
neutron decay products in the ambient plasma, has the 
opposite effect. 



VII. SUMMARY 

In this paper we suggest and analyze a new acceleration 
mechanism, which operates via continuous conversion of 
accelerated particles from charged into neutral state and 
back. The proposed converter mechanism is efficient for 
acceleration of both protons and electrons (positrons). It 
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is capable of producing the highest-energy (> 10^° eV) 
cosmic rays in either GRB or AGN environments. With 
a much lower energy limit, the mechanism can possibly 
operate in microquasars as well. In the regions of high 
optical thickness, the converter mechanism is an efficient 
means of transferring the kinetic energy of bulk relativis- 
tic flow to the accompanying radiation, which could ex- 
plain, for example, the origin of GRB emission. Some 
peculiarities of the accompanying emission can be a tell- 
tale sign of the converter mechanism. For example, the 
production of the highest-energy CRs should be linked 
with the powerful neutrino emission at a level at least 
comparable to the power in CRs. 

Despite a certain similarity to the standard (diffusive) 
acceleration mechanism, the converter mechanism vio- 
lates some of its inherent relations. For example, the 
maximum particle energy attainable in the converter 
mechanism is F^/^ times larger than in the standard one, 
provided acceleration occurs in a shear flow with chaotic 
ambient magnetic field. Also, the maximum energy of 
synchrotron photons appears to be F^ times larger. 



Generally speaking, a beam pattern wider than 1/F 
is characteristic for any type of accompanying emission 
(synchrotron radiation at highest energies, neutrino emis- 
sion or photons from pion decay), as well as for the es- 
caping neutrons. This distinctive feature of the converter 
mechanism opens an interesting possibility for observa- 
tion of the off-axis blazars and GRBs. 
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